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Abstract 



We calculate the free energy of a hot gluon plasma to leading order in hard- 
pg , thermal-loop perturbation theory. Effects associated with screening, gluon 

K^ I quasiparticles, and Landau damping are resummed to all orders. The ultra- 

violet divergences generated by the hard-thermal-loop propagator corrections 
cn ' can be cancelled by a temperature-independent counterterm. The deviation 

^I^ ■ of the hard-thermal-loop free energy from lattice QCD results for T > 2Tc 

Q^ , has the correct sign and roughly the correct magnitude to be accounted for 

^^ I by next-to-leading order corrections. 
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Relativistic heavy-ion collisions will soon allow the experimental study of hadronic mat- 
ter at energy densities that should exceed that required to create a quark-gluon plasma. A 
quantitative understanding of the properties of a quark-gluon plasma is essential in order 
to determine whether it has been created. Because QCD, the gauge theory that describes 
strong interactions, is asymptotically free, its running coupling constant as becomes weak 
at sufficiently high temperatures. This would seem to make the task of understanding the 
high-temperature limit of hadronic matter relatively straightforward, because the problem 
can be attacked using perturbative methods. Unfortunately, the perturbative expansion in 
powers of a^ does not seem to be of any quantitative use even at temperatures that are 
orders of magnitude higher than those achievable in heavy-ion collisions. 

The problem is evident in the free energy JF of the quark-gluon plasma, whose weak- 
coupling expansion has been calculated through order a^J"^ [|1])0- ^^ optimist might hope to 
use perturbative methods at temperatures as low as 0.3 GeV, because the running coupling 
constant as(27rT) at the scale of the lowest Matsubara frequency is about 1/3. However, 
the terms of successive orders in a^"^ in JF form a strictly decreasing series only if as is 
less than about 1/20, which corresponds to a temperature greater than about 10^ GeV. 
At temperatures below 1 GeV, the corrections show no sign of converging, although the 
convergence can be somewhat improved by using Pade approximations [^. It is clear that 
a reorganization of the perturbation series is essential if perturbative calculations are to be 
of any quantitative use at temperatures accessible in heavy-ion collisions. 

The poor convergence of the perturbation series is puzzling, because lattice gauge theory 
calculations indicate that the free energy JF of the quark-gluon plasma can be approximated 
by that of an ideal gas unless the temperature T is very close to the critical temperature 
Tc for the phase transition 0J^ . The deviation of JF from the free energy of an ideal gas of 
massless quarks and gluons is less than about 25% if T is greater than 2Tc. Furthermore, 
the lattice results can be described surprisingly well for all T > Tc by an ideal gas of quark 
and gluon quasiparticles with temperature-dependent masses [0. 

The large perturbative corrections seem to be related to plasma effects, such as the 
screening of interactions and the existence of quasiparticles, which arise from the momentum 
scale a^^T. One possible solution is to use effective-field-theory methods to isolate the effects 
of the scale T 0,^ , and then use a nonperturbative method to calculate the contributions 
from the lower momentum scales of order a^^T and smaller. The effective field theory for 
the lower momentum scales is a gauge theory in three euclidean dimensions, and it can be 
treated nonperturbatively using lattice-gauge-theory methods. Such an approach has been 
used by Kajantie et al. to calculate the Debye mass for thermal QCD p]]. One of the 
limitations of this approach is that it can not be applied to the real-time processes that are 
the most promising signatures for a quark-gluon plasma. 

An analogous convergence problem arises in the free energy of a massless scalar field 
theory with a 0^ interaction, and several approaches to this problem have been proposed 
PJTO[| . One of the most promising approaches is "screened perturbation theory" developed 
by Karsch, Patkos, and Petreczky 0, which involves a selective resummation of higher order 
terms in the perturbative expansion. This approach can be made systematic by using the 
framework of "optimized perturbation theory" [|ri|] . A mass term proportional to 0^ is added 
and subtracted from the lagrangian, with the added term included nonperturbatively and 
the subtracted term treated as a perturbation. The renormalizability of the mass term guar- 



antees that the new ultraviolet divergences generated by the mass term can be systematically 
removed by renormalization. When the free energy is calculated using screened perturba- 
tion theory, the convergence of successive approximations to the free energy is dramatically 
improved. 

A straightforward application of screened perturbation theory to a gauge theory like 
QCD is doomed to failure, because a local mass term for gluons is not gauge invariant. 
However there is a way to incorporate plasma effects, including quasiparticle masses for 
gluons, into perturbation theory in a gauge-invariant way, and that is by using hard-thermal- 
loop (HTL) perturbation theory. This involves adding and subtracting HTL correction 
terms to the action |T^, treating the quadratic parts of the added terms nonperturbatively 



and treating the remaining terms as interactions. The resulting effective propagators and 
vertices are complicated functions of the energies and momenta. The nonlocality of the 
HTL correction terms raises conceptual issues associated with renormalization, since the 
ultraviolet divergences they generate may not have a form that can be cancelled by local 
conterterms. 

In this Letter, we calculate the free energy of a hot gluon plasma explicitly to leading 
order in HTL perturbation theory. In spite of the complexity of the HTL propagators, their 
analytic properties can be used to make calculations tractable. Although complicated ultra- 
violet divergences arise in the calculation, many of them cancel. The remaining divergence 
is removed by a temperature-independent counterterm at the expense of introducing an 
arbitrary renormalization scale. With reasonable choices of the renormalization scales, the 
deviation of the HTL free energy from lattice QCD results for T > 2Tc has the correct sign 
and roughly the correct magnitude to be accounted for by next-to-leading order corrections. 

Our starting point is an expression for the free energy from the one-loop gluon diagram in 
which HTL corrections to the gluon propagator have been resummed. In the imaginary-time 
formalism, the renormalized free energy can be written as 

^HTL = 4(rf - 1)^ \og[ujl + k^ + Ht] + 4^ log[A;2 - H^] + A^, (1) 

where d is the number of spatial dimensions and AjF is a counterterm. The transverse and 
longitudinal HTL self-energy functions are 
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where rrig is the gluon mass parameter. The sum-integrals in (|l|) represent 
TJ2nf^^~'^ I d'^k/i'^'^Yj where the sum is over the Matsubara frequencies Un = 2TmT . If we 
use dimensional regularization to regularize ultraviolet divergences, AJ-" cancels the poles 
in d — 3 in the sum-integrals and /x is the minimal subtraction renormalization scale. If we 
set Ht = n^, = 0, the free energy ([|) reduces to that of an ideal gas of massless gluons: 
^ideai = -(87rV45)T4. 

Standard methods can be used to replace the sums over n in (|l]) by contour integrals in 
the energy uj = iujn- The integrands are weighted by the thermal factor l/{e^^ — 1) and the 



contour encloses the branch cuts on the real cu-axis. The arguments of the logarithms have 
branch cuts associated with Landau damping that extend from —k to +k. The integrands 
also have logarithmic branch cuts that end at the points u = ^uxik) in the transverse term 
and at u = ±ujL{k) in the longitudinal term, where uxik) and oJiik) are the quasiparticle 
dispersion relations for transverse gluons and longitudinal gluons (plasmons), respectively. 
These dispersion relations are the solutions to the following transcendental equations ||13||: 
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By collapsing the contours around the branch cuts, we can separate the integrals over uj into 
quasiparticle contributions and Landau-damping contributions. These individual contribu- 
tions have severe ultraviolet divergences. The divergences can be isolated by subtracting 
expressions from the integrands that render the integrals finite in d = 3 and then evaluating 
the subtracted integrals analytically in d dimensions. If we impose a cutoff A on fc and 
a;, there are power divergences proportional to A^ and m?h? and logarithmic divergences 
proportional to m^T^ log A, m^log^A, and m^logA. The A^ divergence is cancelled by 
the usual renormalization of the vacuum energy density. The m^ log^ A divergences cancel 
between the quasiparticle and Landau-damping contributions to the transverse term. The 
cancellation can be traced to the fact that Ht in (0) is analytic in the energy uo = iuon at 
oj = oo. The temperature-dependent m^T^ log A divergences cancel between the longitu- 
dinal and transverse terms. This cancellation is clearly related to gauge invariance, which 
relates the coefficients of IIt and IIl in (0) and (^. The remaining divergences arise from 
integration over large three-momentum and are cancelled by the counterterm AT in (|l]). In 
dimensional regularization, power divergences are set to zero and the logarithmic divergence 
appears as a pole in d — 3. In the minimal subtraction renormalization prescription, it is 
cancelled by the counterterm AjF = — 9m^/(87r^((i — 3)). 

Our final result for the free energy of the gluon plasma to leading order in HTL pertur- 
bation theory is 
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where /is = v47re^'''/^yU is the renormalization scale associated with the modified minimal 
subtraction {MS) renormalization prescription and 7 is Euler's constant. The first term in 
(^) is the free energy of an ideal gas of transverse gluons with dispersion relation ujt- The 
second term is the free energy of an ideal gas of plasmons with dispersion relation u^, with a 
subtraction that makes it vanish in the high-temperature limit. The third term is a Landau 
damping contribution that involves angles (^l and (^t defined by 
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where K^ = k"^ — uj"^ and L = log[(A; + uj)/{k — u)]. Both 0t and 0l vanish at the upper 
endpoint A; — * 00 of the integral over k. At the lower endpoint k ^ u, (pT vanishes and 
(pL approaches vr. The terms in (^) proportional to m?T'^ and mft come from the zero-point 
energies of the quasiparticles and from subtraction integrals. In the high-temperature limit 



nig <^ T, J^HTL can be expanded in powers of nig/T: 
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where a 
limit T - 



= 3?7i^/(47r^T^). Only integer powers of a appear beyond the a^^^ term. In the 
with nig fixed, Thtl is proportional to m^: 
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This low-temperature limit is sensitive to the value of /is. In particular, the coefficient of 
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in ( ]T0| ) changes sign at /is = 0.717m 
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The free energy of a quark-gluon plasma in the high-temperature limit has been calcu- 
lated in a weak-coupling expansion through order a^J"^ |T]J3]- The result for a pure gluon 
plasma with A^c = 3 is 
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where a = asi^nT)/?!. In the limit a^ — > 0, the gluon mass parameter nig is given by 
m'l = (47r/3)asT^. The expansion parameters a in (^) and (0) therefore coincide in this 
limit. The order-a^^^ terms in these expansions are identical, because HTL resummation 
includes the leading effects associated with Debye screening. Note that HTL resummation 
overincludes the order-a correction by a factor of three. The remaining order-a corrections 
would appear at next-to-leading order in HTL perturbation theory. The order-a correction 
in Thtl together with the corrections that are higher order in a combine to give a total 
correction that is negative, in spite of the large positive contribution from the a^/^ term. 

In this leading-order calculation, T, m^, and /is all appear as independent parameters. 
The parameters nig and /is should be chosen as functions of T and Ug so as to avoid large 
higher order corrections in HTL perturbation theory. At asymptotically large temperatures, 
the fractional correction to Tdeax from the next-to-leading order diagrams must reduce to 
+ (15/2)as/7r in order to agree with Tqcd up to corrections of order a^. This will require 
setting the thermal gluon mass parameter to 

Att 
^%T) = ya,(/i4)T2, (12) 

with a renormalization scale /i4 of order T. A reasonable choice is /i4 = 27rT, the euclidean 
energy of the lowest Matsubara mode. The logarithmic divergences associated with the 



three-dimensional renormalization scale /is will be cut off by higher order corrections at a 
scale of either nig or T. If we choose /is to be of order m^, the a^log(/i3/T) term in (H) will 
reproduce a fraction of the a1 loga^ term in Tqcd- A reasonable choice is /xs = A/3mg(T), 
which is the Debye screening mass. 

In Fig. ^ we compare various approximations to the free energy of a gluon plasma to 
the lattice results for pure-glue QCD from Boyd et al. 0]. The unshaded bands are the 
ranges of the perturbative expansions of the QCD free energy when the renormalization 
scale is varied by a factor of 2 from the central value /i4 = 27rT. The four bands correspond 
to Tqcd in ([TTp truncated after the a^, o?J'^ ■, cP'si ^^^d ot'J'^ terms, respectively. We use 
a running coupling constant that runs according to the two-loop beta function: as(/i4) = 
(47r)/(lll)[l- (102/121) log(L)/L], where L = \og{fj,l/ Aj^) . The parameter Ajg^ is related 
to the critical temperature T^ by Tc = LOSA^g^ [|]. The poor convergence properties of the 
perturbative expansion and the strong dependence on the renormalization scale are evident 
in Fig. |I]. The shaded region in Fig. |I] is the range of the HTL free energy Thtl when 
the renormalization scales are varied by a factor of 2 from the central values /i4 = 2'kT and 
/is = V^rrigiT). For these choices of /is and /i4, Thtl/T'^ is a slowly increasing function of T. 
This feature follows from the fact that mg{T) is approximately linear in T, with deviations 
from linearity coming only from the running of the coupling constant. If mg{T) was exactly 
linear in T, Thtl/T'^ would be independent of T. With our choices of /is and /i4, the 
HTL free energy lies significantly below the lattice results for T > 2Tc. This should not be 
of great concern, because the next-to-leading order correction in HTL perturbation theory 
will give a fractional correction to J-'ideai that approaches +(15/2)a;s(/i4)/7r at asymptotic 
temperatures. It has the correct sign and roughly the correct magnitude to decrease the 
discrepancy with the lattice QCD results at the highest values of T. With the inclusion of the 
next-to- leading order correction, the error at asymptotic temperatures will fall like al log as- 
If the next-to-next-to-leading order correction was also included, the error would decrease 
to order a^ logofs. Because of the magnetic mass problem, the error can be decreased below 
order a^ only by using nonperturbative methods. 

We have proposed HTL perturbation theory as a resummation prescription for the large 
perturbative corrections associated with screening, quasiparticles, and Landau damping. 
The free theory around which we are perturbing is similar to the phenomenological quasipar- 
ticle models, but the effects of interactions between the quasiparticles can be systematically 
calculated. This approach can be applied to the real-time processes that may serve as sig- 
natures for the quark-gluon plasma. We have demonstrated that HTL perturbation theory 
is tractable by calculating the leading term in the free energy of a pure gluon plasma. With 
reasonable choices of the renormalization scales, the deviation of the hard-thermal-loop free 
energy from lattice QCD results for T > 2Tc has the correct sign and roughly the correct 
magnitude to be accounted for by next-to-leading order corrections. A challenging problem 
is to extend the calculation of the free energy to next-to-leading order in HTL perturbation 
theory. If the next-to-leading order correction proves to be small for temperatures within 
an order of magnitude of T^, we may finally have a perturbative framework that will al- 
low quantitative calculations of the properties of a quark-gluon plasma at experimentally 
accessible temperatures. 
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FIG. 1. The free energy of a hot gluon gas norniahzed to that of an ideal gas of gluons as a 
function of T/T^- The black diamonds are the lattice results for pure-glue QCD from Ref. Q. The 
bands enclosed by the curves labelled 2, 3, 4, and 5 are the QCD free energy Tqcd truncated 
after the a^, aj , a^, and a J terms, respectively. The bands correspond to varying ^4 by a 
factor of 2 from the central value 27rT. The shaded region is the HTL free energy J-htl with 
mi = (47r/3)Q;s(M4)^^- The region corresponds to varying the renormalization scales by a factor 
of 2 from the central values //3 = ^/^mg and ^4, = 2t:T. 



